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Abstract. The deflection angle A(j) of light rays by the 
gravitational field of a spherical system M(r) is calculated 
using the Modified Newtonian Dynamics (MOND). It is 
shown that A(f> with an impact parameter tq can be ex- 
pressed by the measured rotation velocity v{r) as 



A0(ro) 



where 



v{r] 



'{Ga,M{r)f'\ 

' GM{r) \ 



and Tc is the critical radius that is determined by the 
critical acceleration aq. In the Newtonian limit of the 
gravitational acceleration a ^ ao, approaches A(j) = 
2Gm{ro) / c^ro with the projected surface mass m(ro). 
Whilst the asymptotic value of Acj) reaches a constant 
7r(wcx)/c)2 in the low-acceleration limit of a ^ oq. Taking 
the empirical correction of a factor of 2 from the theory 
of general relativity into account and utilizing the rela- 
tion between rotation velocity v and velocity dispersion 
(7, MOND results naturally in a constant defiection an- 
gle, 47r(a/c)2, which has been widely used in the present- 
day study of gravitational lensing by galaxies and clusters 
of galaxies, implying that without introducing the mas- 
sive halos acting as r"^ for dark matter MOND has no 
difficulty in reproducing the known cases of gravitational 
lensing associated with galaxies and clusters of galaxies. 
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1. Introduction 

Many attempts have been made for decades to account 
for the large discrepancies between the luminous matter 
seen from the galaxy radiations and the dynamical matter 
derived from galaxy rotation curves. A conventional way 
out of these mass discrepancies is to assume the existence 
of some kiirds of unseen matter, which has been widely ac- 
cepted today by the community of astronomers and physi- 
cists and many experiments have been devoted to searches 
for these dark matter despite that there have been no any 
positive results obtained thus far. Whilst some alterna- 
tives to the puzzle are to modify the Newtonian inverse 
square law, among which the most successful of the empir- 
ically suggested models is the Milgrom's MOdified New- 
tonian Dynamics (MOND) (Milgrom 1983a, b, c). 

As a proposal to avoid introducing dark matter in as- 
trophysics, MOND has naturally predicted the flat form 
of the rotation curves for spiral galaxies, the TuUy-Fisher 
law and the luminosity-rotation velocity/velocity disper- 
sion relations in spiral/elliptical galaxies. Although there 
have been some arguments, besides the unknown details of 
any fundamental physics, against MOND using the stel- 
lar kinematics in dwarf galaxies (Lake 1989; Gerhard & 
Spergel 1992) and the dynamical properties of X-ray se- 
lected clusters of galaxies (Gerbal et al. 1992), these claims 
still need to be further considered (Milgrom 1991; 1993). 

As was pointed out recently by Sanders and Begeman 
(1994), the more fundamental problem with MOND is the 
observation of gravitationally distorted images of back- 
ground galaxies by the gravitational potential of interven- 
ing clusters of galaxies. Indeed, the deflection of light rays 
by the gravitational field has remained to be an unsolved 
problem in MOND (Milgrom 1993, private communica- 
tion). The presently known gravitational cases, such as 
the lensed quasar pairs, luminous arcs (see Surdej & Sou- 
cail 1993 for a summary), quasar-galaxy associations (see 
Wu 1994 for a summary), etc., associated with galaxies 
and clusters of galaxies where MOND would appear to be 
important can constitute a critical test for MOND. 
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This paper is then to calculate the deflection of light 
by the gravitational field of a spherical matter system in 
MOND. Due to the absence of a general theory of gravity 
in MOND comparable to General Relativity, the present 
paper deals with the motion of particles only in its classical 
way. 

2. The equation of motion in MOND 

In MOND the true gravitational acceleration g is related 
to the Newtonian gravitational acceleration gn as 



M(5/ao)g = gn, 



(1) 



where oq is a new fundamental constant or the critical 
acceleration parameter with the value of oq ~ 10~^ cm 
s~^ (Milgrom & Braun 1988), and fj, is some function that 
has the asymptotic behaviour 

(i{x) = 1, a; » 1; 
IJ.{x) = X, X <^1. 



In particular, the major properties of the results in MOND 
are insensitive to the choice of /x (Milgrom 1983b). For 
simplicity, we take fi{x) = 1 when x > 1 and fi{x) = x 
when a; < 1. 

The motion of a test particle in the gravitational field 
of spherical symmetry is described by 



.«0M(f:) = 0. 



GM(r) 




(3) 



Here and represent the components of acceleration 
along the r- and (f)- direction, respectively. M(r) is the 
total mass enclosed within the radius of r. The equation of 
motion of the test particle can then be specifically written 
as 



r — r6'^ = —^ 



r'^0 = ci; 



GaoM{r) 



r > Tc 



(4) 



and 



r'e = C2; 



r < Tr. 



(5) 



where fc measures the radius separating Newtonian me- 
chanics from MOND and is determined by 



Igm(tc 

ao 



(6) 



two different kinds of equations of motion reflects their 
significant difference. Setting r = f = yields 



v{r) 



\GaoM{r)f\ r > r,; 

1/2 



(7) 



r < rr. 



For a pointlike mass or a finite mass distribution within 
a radius of R, a constant circular velocity instead of the 
decreasing form of r~^/^ in the Newtonian mechanics is 
predicted by MOND when r tends to infinity. This pro- 
vides a natural explanation for the flatness of the rotation 
curves of spiral galaxies alternative to dark matter. 

Some initial and/or boundary conditions are needed to 
find the solution to the orbit of a test particle from eqs.(4) 
and (5). Define an impact parameter ro 



(2) r(ro) = 



so that the velocity of the test particle at r 
written as 

vo = ro9{ro) 



(8) 



ro can be 



(9) 



Moreover, the first-order time derivatives of position of the 
test particle at the critical radius Tc are required to keep 
their continuities when the particle crosses the boundary 
oi r = Tc from either sides. An immediate consequence of 
this requirement is that the two integral constants Ci and 
C2 in eqs. (4) and (5) have the same value and 



ci = C2 = voro 



(10) 



In fact, this is the result of the conservation of angular 
momentum analogue to the one in Newtonian dynamics. 

The equation of motion of the test particle, eqs. (4) and 
(5), can be unified to one simple equation by replacing the 
matter distribution M(r) by the circular rotation velocity 
v{r) of eq.(7) 



1-1^ 



v{r) 



= 0. 



(11) 



The solution to this equation with the boundary condition 
of eq.(9) is 



dr 



= ±- 



1 



V^' - r'o ^ 



1 



Jr. 



2t;(r)2 
,5 



dr 



(12) 



Ci and C2 are the integral constants, and the time deriva- 
tive of the position of the test particle is denoted by 
The circular rotation velocity v(r) resulting from these 



Given the observed rotation curve v(r) in and around the 
massive system or the matter distribution M(r), the orbit 
of a test particle can be found from eq.(12). 
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3. The deflection of light in MOND 



The deflection angle of light rays propagating from and to 
the distant universe with an impact parameter of ro from 

the massive body M is 



A(/) = 2|^;.(ro)-^(oo)| -TT. 



(13) 



The test particle discussed in the above section is now 
replaced by the photons and hence, the speed of light c is 
used instead of vq. 

Eq.(12) can now be expanded in the series of {v{r)/cf: 



<P{r) 



± 



rodr 



(r2 - r^)3/2 



0{ 



v{r) 



(14) 



'(r) dr 
72 



approaches infinity 



Note that the integral of 

^ Jrn r ^ ^ ^ 

when r — > oo, which may then invalidate the expansion 
of cq.(12) into eq.(14). To see this, utilizing a constant cir- 
cular rotation velocity v^o at a large distance (r*) well be- 
yond both the radius of the massive body and the MOND 
limit leads to 



2 dr 

r 



c / r* 



However, the distance r that satisfies r ~ r* exp{{c/voo)^) 
would exceed the whole visible size of the present universe 
for any realistic gravitationally bound systems like galax- 
ies and clusters of galaxies. Therefore, eq.(14) is valid in 
the presently observable universe. The first term of the 
right-hand side of eq.(14) is the orbit of a straight line 
when the photon travels in an Euclidean space, and the 
rest terms represent the contributions of the gravitational 
field. 

Finally, to the first order of {v{r)/cY the deflection of 
light is 



'■f 

Jrn 



ro 



ry/r'^ — ■ 



zdr. 



(15) 



3.1. Pointlike 



In the case of pointlike mass, inserting eq.(7) with M(r) = 
M into eq.(15) yields: 
(l)For ro > rc. 



VGaoM 



(16) 



(2)For ro < r, 

2GM 



c^ro 



ro 2GM 



rc + ro 



ro 



rc + ro 



, 2^/GaoM . _i ro 
H 7i sm — . 



In particular, eq.(17) reads 



2GM 
c^ro 



00, 



(17) 



(18) 



i.e., the result of Newtonian limit. Whilst taking rc — ro 
in eq.(17) identifies the result in a non-Newtonian limit 
[eq.(16)]. 

The most interesting result appears in the low- 
acceleration limit or at large radii of galaxies and clusters 
of galaxies where MOND plays an important role. Eq.(16) 
can be simplified using the observable asymptotic rotation 
velocity Wqo 



(19) 



The significant feature of the defiection of light in MOND 
is its constant angle of light bending, in accord with the 
flat rotation curves at large distances from the centers 
of galaxies and of clusters of galaxies. A relatively larger 
deflection is then provided by MOND than that by the 
Newtonian mechanics (see Fig.l). 



Fig. 1. The deflecting angle of light A0 in unit of 
(po = 27r(woo/c)^ by a pointlike mass versus the impact param- 
eter ro in unit of Vc. Beyond the critical distance r^, MOND 
predicts a constant deflection of light instead of the declining 
form of 1/r derived from the Newtonian dynamics (the dotted 
line) 



Recall that the present theory of gravitational lens- 
ing uses indeed a constant deflection of light for the dark 
matter associated with galactic halos (Turner, Ostriker 
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& Gott, 1984) and for matter distributions in clusters 

of galaxies (Wu & Hammer, 1993; Wu 1993). Essentially, 
these constant deflections of light can account for the sta- 
tistical properties of the lensed quasar pairs and the gi- 
ant luminous arcs. To quantitatively compare the result 
of MOND with the update theory of gravitational lensing 
that adopts the line-of-sight velocity dispersion for galax- 
ies and clusters of galaxies, the circular velocity v^c in 
eq.(19) should be written in terms of the observed line- 
of-sight velocity dispersion. Unfortunately, there has not 
been set up an analog of the virial theorem in MOND (Mil- 
grom 1983c). Additionally, the present discussion is con- 
fined to a pointlike mass rather than an extended system, 
although the point mass can be regarded as a good ap- 
proximation to a galaxy seen at a large distance. Note that 
there is no massive dark halo at all surrounding the lumi- 
nous galaxy according to MOND. Basically, in a singular 
isothermal sphere the circular velocity is related to the 
line-of-sight velocity dispersion by v = \/2cr. Therefore, 
the gravitational deflecting angle by a singular isothermal 
sphere used in gravitational lensing today can be written 
as 



477^ = 277^. 



(20) 



Here both v and cr in a singular isothermal sphere are 
independent of radius and thereby, v = Voc ■ This property 
is the unique among the various matter distributions (Wu 
1993). 

Due to the absence of the general theory of gravita- 
tion in MOND, a "quasi-general relativity" correction of 
a factor of 2 to the light bending derived from a classical 
way in MOND is empirically adopted. As a consequence 
of this, eq.(19) reads 



specified: 
(l)For ro>rh, 



A(j) = TT- 



.2 ' 



(2)For ro < rn, 



A(j) = 2- 



sm 



-1 ro 



+ 



J To \ ^OO 



rodr 



(22) 



(23) 



Figure 2 illustrates the variations of with the impact 
parameters for three kinds of rotation curves. Note that 
utilizing the observed rotation curves in the calculation of 
is independent of the assumed dynamics of whether it 
is the Newtonian one or the MOND. 



A</) = 27r^, 



(21) 



which reconcile the presently used formula in the study of 
gravitational lensing by larger systems (galaxies, clusters 
of galaxies, etc.). The key point is that MOND provides a 
deflection of light identical to the one by a singular isother- 
mal sphere without introducing the massive halos having 
the density form of r~^. It then turns out that all the 
presently known cases of gravitational lensing associated 
with galaxies and clusters of galaxies can be explained by 
MOND using the hmiinous matter alone that concentrated 
in the central regions. 

3.2. Extended mass distribution: rotation velocity 

For a galaxy whose circular velocity has been well mea- 
sured, A(j) can be straightforward obtained using eq.(15). 
A typical rotation curve exhibits an increasing form in the 
central region and tends to a constant Voo beyond a char- 
acteristic radius of rh- Two types of light orbit should be 



Fig. 2. The deflecting angle of light A0 in unit of 
(po = 7r(voo/c)^ derived from the observed circular veloc- 
ity. Three types of rotation curves in the central region 
of radius of are assumed: (l)w(r) = Voc{r /thY^^ (top); 
(2)t)(r) = Vaair/rn) (middle); (3)ti(r) = Voo{r/rhY (bottom). 
Beyond the radius of rh the rotation velocity remains un- 
changed, leading to a constant deflection of light. 



3.3. Extended mass distribution: surface brightness 

Surface brightness (optical/X-ray) of galaxies and of clus- 
ters of galaxies are relatively easier to measure than their 
rotation velocities. Assuming that the light profile traces 
the mass distribution as argued By MOND, one can ob- 
tain the deflection of light from eqs.(7) and (15): 
(l)For xq < Xc, 



A(l){xo) 



XcXa 



\/m{xc) 



Xo 



M{x)dx 

"^sj X^ — X 
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(2)For Xo > Xc 



^M{x)dx 



(24) 



A</)(a:;o) 



■I xn 



M(x)dx 



Xl/x^ — : 



(25) 



Here (^o = ^^fGaoM^jf?, and M(x) is the dimensionless 
mass enclosed within x normalized to the total mass Mt. 
All the distances [x) are measured in unit of the so-called 
length scale in the model fit to the surface luminosity dis- 
tribution. The subscripts "c" and "o" indicate the critical 
and the impact positions, respectively. 

An example is shown in Fig. 3 using the King model 
for the matter distribution. This model has been found to 
fit fairly well the surface brightness of the central parts of 
galaxies and of clusters of galaxies, and has been widely 
adopted today in the modeling of the luminosity distribu- 
tion for these systems. The spatial mass distribution M(r) 
can be found to be 

M(x) = ^(ln(x + yrT^)--^), 

in which x = r/b and h is the core radius. The deflection 
of light resulting from this model shows a slight increasing 
form instead of the declining one of ~ in the Newto- 
nian limit, when the impact distance tends to infinity. This 
shows again that MOND provides a significantly large de- 
flection of light beyond oq.. 



4. Discussion and conclusions 

Dark massive halos of galaxies and unseen matter excess 
in clusters of galaxies, if they were real, can be well de- 
scribed by a singular isothermal sphere with a power law 
of r~^. This model has been widely used in the study of 
gravitational Icnsing today. Certainly, this is due partially 
to the fact that it results in a constant deflection of light 
47r(c7/c)^, making the calculations much simplified. Nev- 
ertheless, there have been no any conclusive evidences, 
from the experiments that search for dark matter, for the 
existence of the hidden mass. Dynamical analysis of the 
rotational curves of galaxies on the basis of Newtonian 
mechanics is the only strong observational evidence that 
shows the mass discrepancies. 

An alternative to the missing mass is to modify the 
Newtonian gravity. Though it is not conventional, Mil- 
grom's MOND has successfully predicted the flat rotation 
curves of galaxies, leading to a new sight into the dynamics 
on large scales. The present paper has computed one of the 
critical issues in MOND, the deflection of light rays by a 
spherical gravitational field. Surprisingly, MOND provides 



Fig. 3. The deflection of light A(j)/<j)o by the matter distribu- 
tion following a King model. Four critical radii in unit of core 
radius b have been shown: (from top to bottom) rc = 1006, 
rc — 10b, rc = 56 and rc = 6. The normalization <f>Q is 
(9aoba^y^^/c^. The results from the Newtonian dynamics are 
also plotted in the dashed lines for comparisons. 



a constant deflecting angle at large distance from the cen- 
ter of the gravitational field, which is consistent with the 
the value presently used in the study of gravitational lens- 
ing by galaxies and clusters of galaxies. It is then likely 
that all the lensing cases can be equally reproduced in 
MOND without the massive dark matter in galaxies and 
in clusters of galaxies. Similar to the flatness of rotation 
velocity in galaxies predicted by MOND without assum- 
ing the massive halos, the constant deflection of light from 
MOND has the same effect as the halos did. 

It has been shown that light bending is no more a crit- 
ical argument against MOND today. Conversely, MOND 
predicts a reasonable deflection angle of light by large mas- 
sive systems. Therefore, whether or not MOND reflects 
the nature of gravity needs to be further investigated us- 
ing other astronomical methods. 
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